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Abstract.  To  estimate  the  genetic  component  of  cellular 
membrane  thermostability  in  Phaseolus  vulgaris,  parental 
and  Fj^  plants  from  a  6-parent  half  diallel  cross  were 
tested  by  electrical  conductivity  and  the  results  were 
analyzed  by  Hayman's  method.  Membrane  thermostability 
was  found  to  be  a  quantitative  trait,  with  environmental 
and  dominance  effects  accounting  for  most  of  the 
phenotypic  expression.  Narrow  sense  heritability  was 
low  (5.4%). 
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It  is  well  known  that  high  temperature  is  one  of 
the  major  factors  limiting  yield  of  common  bean.  In 
order  to  breed  for  heat  tolerance,  one  needs  information 
about  the  inheritance  of  the  character  and  a  suitable 
method  of   screening   for  it. 

Screening  methods  and  inheritance  of  heat 
tolerance  in  common  beans  have  been  studied  by  several 
researchers.  Benepal  and  Rangappa  (3)  screened  5380 
accessions  for  their  ability  to  set  pods  in  the  field. 
Wien  and  Hunger  (31)  also  tested  lines  for  heat 
tolerance  in  the  field.  Weaver  et  al.  (29)  suggested  a 
possible  screening  procedure  testing  pollen  growth  in 
vi.t_ro.  or  i.n  vivo.  Ng  and  Bouwkamp  (21)  rated  more  than 
600  accessions  according  to  their  response  to  high 
temperature  in  the  greenhouse.  Bouwkamp  and  Summers  (4) 
reported  on  the  inheritance  of  combined  temperature  and 
drought  stress  resistance  based  on  the  number  of  pods 
set   per  plant. 

In  addition  to  morphological  characters,  methods 
of  measuring  cell  viability  have  been  tested  in  the  hope 
of  finding  a  rapid  laboratory  screening  procedure.  The 
electrical  conductivity  method  is  the  most  important  of 
these.  It  tests  for  cell  membrane  thermostability  as 
measured   by   leachate   conductivity   and   has    recently  been 
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used  as  an  index  of  viability  after  heat  treatment. 
Various  plant  species  and  parts  have  been  examined  by 
the  conductivity  test.  Laminar  pieces  of  tomato  and 
epidermal  strips  from  onion  bulb  were  tested  by  Onwueme 
(22),  holly  root  cells  by  Ingram  (12),  and  pear 
suspension  cells  by  Wu  (32).  A  diversity  of  turf 
grasses  has  been  tested  by  with  this  method  (28,30);  and 
Chen  et  al.  (6)  have  examined  tomato,  soybean,  and 
potato. 

Comparison  of  the  electrical  conductivity  test 
with  other  methods  and  with  field  performance  suggests 
that  the  electrical  conductivity  method  is  a  moderately 
reliable  and  convenient  screening  procedure  for 
measuring  heat  tolerance.  Marsh  et  al.  (18)  measured 
heat  tolerance  in  common  bean  by  the  conductivity  test, 
a  hot  water  dip,  percent  pod  set,  and  pollen 
8 t a ina b i 1 i t y.  They  concluded  that  the  conductivity 
method  best  combined  reliability  with  early  testing 
convenience.  Marsh  et  al.  (17,19)  reported  for  the 
conductivity  method  that  the  killing  times  for  the  5 
genotypes  used  were  in  agreement  with  previous  ratings 
from  yield  data.  They  also  studied  heat  tolerance 
inheritance  by  conductivity  tests  for  the  parents,  F^, 
backcross,    and   F2   populations   from  3   crosses   of  heat- 


tolerant  X  hea t- int o 1 e r an t  lines.  Schaff  (23)  observed 
significant  correlation  of  electrical  conductivity  with 
field  performance  under  heat  stress,  developed  a 
sigmoidal  model  to  calculate  killing  temperature,  and 
conducted  a  6-parent  weighted  diallel  analysis  to 
determine  the  inheritance  of  heat  tolerance.  In  soybean 
and  sorghum,  good  correlation  has  been  found  between 
heat  tolerance  as  measured  by  electrical  conductivity 
and   field   performance   measured   by  yield  (20,  25). 

The  Hayman-Jinks  diallel  cross  has  been  widely 
used  for  inheritance  studies  (5,  14),  despite  the  fact 
that  some  reports  have  been  openly  critical  of  the 
diallel  analysis  as  a  method  for  studying  the  genetics 
of  complex  traits  or  as  a  tool  in  plant  breeding  (8). 
Johnson  (14)  pointed  out  two  major  advantages  that  the 
diallel  cross  provides:  compared  to  other  methods 
available,  the  diallel  cross  technique  permits  a  more 
systematic  approach  to  large  scale  studies  of  continuous 
variation  and  a  better  diciplined  analysis  of  the 
resulting  data;  and  the  overall  analysis  provides 
reliable  genetic  information  on  dominance  and 
r e c e s s i V e n e s s  and  on  complementary  non-allelic 
interaction. 

The   diallel    analysis    makes   it   possible    to  predict 
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the  phenotypes  of  the  completely  dominant  parents,  which 
in  turn  suggests  the  possible  limit  of  selection  among 
genes  showing  dominance.  Baker  (2)  emphasized  the  fact 
that  similar  information  could  be  obtained  from 
different  methods  of  analyzing  diallel  crosses,  such  as 
those  developed  by  Griffing  (9)  and  Gardner  and  Eberhart 
(7).  The  assumptions  required  for  the  genetic 
interpretation  in  se 1 f -po 1 1 ina t ing  plants  were  evaluated 
by  Sokol  and  Baker  (24).  Jones  (15)  modified  the 
Hayman-Jinks  method  so  that  it  can  be  conducted  without 
reciprocal   crosses    (half-dial  lei  analysis). 

The  objective  of  this  experiment  was  to  carry  out 
a  six  parent  half-diallel  analysis  to  determine  the 
inheritance   of   heat  killing   time   in  common  bean. 
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Material  and  Methods 

Three  of  the  parents  used  in  this  study  are 
cultivars  previously  reported  to  be  heat  tolerant  (23, 
31):  PI  3  24607  (P2),  ND  364  (P4),  Wyoming  166  (P5).  The 
other  three  are  heat  intolerant  (23):  PI  271998  (P^), 
Oregon  1604  (P3),  Valley  (Pg).  Valley,  Wyoming  166,  and 
ND  364  were  obtained  originally  from  M.  LeBaron, 
University  of  Idaho,  Kimberly;  PI  271998  and  PI  324607 
from  the  USDA  Plant  Introduction  Station,  Pullman,  WA; 
and  Oregon  1604  from  the  Idaho  Seed  Bean  Co.,  Twin 
Falls,  ID.  All  cultivars  have  been  maintained  by  single 
seed  descent   for   at   least   four  generations. 

The  half-diallel  cross  was  made  in  winter  1984. 
The  seeds  of  each  parent  were  sown  in  5.5  X  5.5  X  5.0  cm 
pots  containing  a  potting  mixture  of  vermiculite,  peat, 
perlite,  and  soil;  and  the  pots  were  put  in  a  growth 
chamber  set  for  30  C  and  a  16-hr  light  period  to  ensure 
rapid  and  uniform  germination.  Ten  days  from  seeding, 
the  seedlings  were  transplanted  into  one-gallon  pots 
containing  the  same  soil  mixture  as  in  the  seeding  pots. 
The  plants  were  maintained  in  a  greenhouse  set  for  28/22 
C  day/night  temperature  and  supplementally  lighted  to 
approximate  a  12-hr  photoperiod  to  ensure  uniform  timing 
of   flower   initiation.      The  crosses  were  made  by  hand 
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pollination  as  soon  as  the  plants  began  flowering. 
Cultivars  with  dominant  marker  genes  were  used  as  male 
parents   whenever  possible. 

In  fall  and  winter  of  1  985,  4  runs  of  the 
conductivity  test  were  carried  out,  each  containing  1 
plant  for  each  of  the  21  accessions  (15  Fj^  hybrids  and  6 
parents).  In  each  run,  2  to  3  seeds  of  each  genotype 
were  sown,  transplanted,  and  maintained  in  the 
greenhouse  as  previously  described.  One  plant  of  each 
genotype  was  tested  for  membrane  stability  in  a  random 
sequence.  The  time  schedule  for  each  run  is  listed  in 
Table  1. 

The  plants  to  be  tested  were  acclimated  at 
flowering  stage  for  24  hours  in  growth  chambers  set  for 
a  constant  37. 5C  and  a  16-hr  photoperiod  at  900  M E  sec"^ 
m  .  Immediately  after  acclimation  the  young,  fully 
expanded  leaves  were  picked  for  testing  by  the  procedure 
of  Kinbacher  (16)  with   the   following  modifications. 

Leaf  discs,  one  cm  in  diameter,  were  washed  with 
deionized-dist il led  water,  changed  3  times,  and  put  into 
test  tubes,  5  per  tube,  each  containing  1  ml  of  water. 
The  treatment  tubes  were  kept  in  a  water  bath  set  at 
47C  for  30,  60,  90,  1  20,  1  50  ,  and  180  minutes,  for  each 
of   the  accessions  with  3   replications   for   each  time. 


Control  discs  were  held  at  room  temperature.  After  the 
tubes  cooled,  20  ml  of  d e i o n i z e d - d i s t i 1 1 e d  water  were 
added  to  each  and  they  were  incubated  at  10  C  for  24 
hours.  The  first  conductivity  was  determined  at  25C 
after  incubation.  All  of  the  tubes  were  put  in  boiling 
water  for  15  minutes  to  kill  the  cells  completely. 
After  24  hours  at  room  temperature,  the  second 
conductivity   reading  was   taken  at  2  5C. 

Relative  leakage,  or  injury,  was  calculated  using 
the  equation: 

Relative    leakage    =    1    -   [  ( l-(     /  C2 ) )  /  ( l-(       /       )  )  J  . 
where  Cj^   =  treatment   first   conductance,    C^^   =  control 
first   conductance,    C2  =  treatment  second  conductance, 
and  =  control   second  conductance. 

The  killing  time  (Time50)  is  defined  as  the  time 
at  which  50%  cell  were  injured.  Viability  was  estimated 
using   the   sigmoidal  equation: 

Viability  =     1   /    [1   +  e'^^^i^e-^imeSO)   ]    +  e 

where  B  is  a  rate  parameter  and  e  is  the  deviation  from 
the   regression    line  (27). 

The  error  term  for  the  diallel  analysis  was  the  MS 
ERROR,  the  mean  square  interaction  of  cross  X  run,  which 
can    be    obtained    by    fitting    the    sigmoidal    model  to 
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various    sets   of    data    (Table  2). 

The  diallel  weighting  and  analysis  procedure 
followed  that  of  Schaff  (23)  with  modifications  and  is 
given  in  the  appendix.  The  same  notation  was  used  as  by 
Hayman  (10).  The  caculations  were  conducted  using  the 
SAS  computer  language. 


11 

Results   and  Discussion 

The  Hayman-Jinks  diallel  model  assumes:  parental 
homozygosity,  diploid  segregation,  no  differences 
between  reciprocal  crosses,  no  multiple  alleles,  no 
epistasis,    and  no  linkage  between  the  genes  studied. 

The  first  three  assumptions  were  confirmed  to  be 
valid  by  observation  of  parental  and  hybrid  phenotypes. 
The  remaining  three  were  tested  with  the  methods  of 
Hayman  (10)  and  Jinks  (13),  namely,  uniformity  of  - 
Vj..       The    t-test    for   heterogeneity   of  -   Vj.   is  not 

significant  (  t=0.046,  P=0.97  ),  and  Figure  1  shows  that 
the  regression  slope  of  W^.  on  V^.  is  not  significantly 
different  from  unity  (for  b=l,  P  =  0.30).  Consequently, 
there  is  insufficient  evidence  to  say  that  Hayman's 
model   is   inappropriate   for   this  experiment. 

The  mean  killing  time  (Time50)  and  its  coefficient 
of  the  standard  error  (K)  for  each  genotype  are 
presented  in  Table  3.  Due  to  the  large  standard  errors 
of  TimeSO  and  the  small  mean  square  residuals  of  the 
sigmoidal  model,  K  values  are  relatively  large.  The 
relative  position  of  the  parents  for  heat  tolerance  is 
similar    to    Schaff's    results  (23). 

Since  only  six  parents   were   intercrossed   in  the 
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present  study  (less  than  ten  parents),  the  genetic 
components  estimated  are  appropriate  only  for  this 
particular  set  of  parents,  rather  than  for  the  entire 
population  (10,11). 

The  main  components  of  genetic  variance  are  listed 
in  Table  4,  and  the  estimates  of  heritability  parameters 
are  given  in  Table  5.  Since  (Hj^/D)^/-^,  which  estimates 
the  degree  of  dominance,  is  larger  than  1,  overdominance 
exists,  as  also  indicated  by  the  negative  intercept  in 
Figure  1.  Also,  the  correlation  coefficent  between  the 
parental  order  of  dominance  (W^+V^)  and  the  weighted 
parental  values  is  very  close  to  0  (r=0. 00008,  P=0.99), 
indicating  that  there  are  equal  numbers  of  positive  and 
negative  genes  showing  dominance.  Marsh  et  al.  (19)  also 
found    that    the  mean    killing    time    exceeded  the 

midparent  for  all  3  crosses  studied  and  interpreted  the 
fact  as  gene  interaction.  Since  their  data  also  fit  the 
additive-dominant  model  with  small  epistatic  effects, 
this   "gene    interaction"  means  dominance. 

For  our  results,  the  ratio  of  dominant  to 
recessive  genes  equals  1.40,  indicating  that  there  are 
more  dominant  than  recessive  genes  for  heat  tolerance. 
There  is  at  least  one  gene  group  showing  some  degree  of 
dominance,    as   indicated   by  the  estimate,    h^/H,    =  0.52. 
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The  proportion  of  genes  with  positive  and  negative 
effect  in  the  parents,  Hj/AH^,  is  0.23,  possibly 
indicating  slightly  unequal  distribution  of  positive  or 
negative   genes   among   parents.  '5 

Significant  differences  were  found  among  runs 
(Table  6),  which  suggests  that  unknown  environmental 
factors  affected  plants  by  causing  different  responses 
to  the  heat  stress  among  runs  even  for  the  same  cross. 
Tal  and  Shannon  (26),  also  using  the  leaf  disc 
conductivity  test  for  membrane  heat  stability,  found 
that  all  Lycopersicon  and  So  1 anum  species  tested  had 
more  injury  in  the  winter  than  in  the  summer.  The 
results  of  our  experiment  support  this  because,  as  shown 
in  Table  1,  overall  means  of  each  run  decreased  as  the 
treatment  date  shifted.  So  weather  conditions,  such  as 
light  intensity  and  daylength,  may  account  for  some  of 
the  variation  among  runs,  among  crosses,  and  among 
plants  of  the  same  cross.  In  order  to  obtain  a  more 
accurate  estimate  of  the  inheritance,  the  present 
conductivity  procedure  needs  to  be  refined  to  minimize 
environmental  error.  Additional  replication  of 
genotypes  tested,  both  within  and  between  runs,  would 
give  more   accurate   estimation  of  killing  times. 
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The  estimate  of  narrow  sense  heritability  is  low 
(5.4%),  which  contrasts  with  the  high  heritability  (59%) 
observed  by  Schaff  (23).  Several  factors  should  be 
taken  into  consideration  when  comparing  these  two 
results.  Firstly,  some  of  the  bean  lines  he  used 
differed  from  those  of  this  study.  Secondly,  his 
experiment  was  an  unbalanced  split  plot  design  using 
killing  temperature,  not  killing  time.  Thirdly,  his 
original  data  were  not  homogeneous  for  Wj.-V^  until  one 
parent,    Oregon  1604,    was  removed. 

By  comparing  the  estimate  of  broad  sense 
heritability  (34.6%)  and  that  of  narrow  sense 
heritability,  it  is  clear  that  the  dominant  effects 
accounted   for   most   of   the  genetic  variation. 

The  narrow  and  broad  sense  her itab il it ie s  of  heat 
tolerance  calculated  by  Marsh  et  al.  (19)  for  their 
different  parental  crosses  ranged  from  2.9%  to  24%  and 
0.0%  to  21.6%,  respectively.  Our  narrow  sense 
heritability  (5.4%)  is  within  their  range,  and  its  low 
value  is  due  to  the  large  environmental  and  dominant 
effects  and  the  small  additive  effects.  Also,  the 
estimate  of  the  heritability  in  this  experiment  was 
based  on  individual  plant  responses,  where  large  errors 
of   estimation  were  common. 
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We  conclude  that  membrane  thermostability  in 
common  bean  is  a  quantitative  trait  and  easily 
influenced  by  the  environment.  Dominant  effects 
accounted  for  most  of  the  genetic  variation.  Because  of 
the  low  heritability  of  cellular  membrane 
thermostability,  should  the  conductivity  method  be 
employed  in  a  breeding  program,  continuous  evaluation 
and  selection  using  large  samples  will  be  required  in 
later  generations.  Even  then,  the  time-  and  labor- 
consuming  nature  of  this  procedure  makes  it  difficult  to 
use.  Unless  the  test  can  be  improved  in  these  respects, 
alternative  testing  procedures  will  have  to  be  developed 
for  an  effective  heat-tolerance   breeding  program. 
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Table  1.  Time  schedule  and  overall  mean  Time50  for  each 
run. 


Run  Planting  Testing  Mean 

1  09/25/85  11/06/85  120.4 

2  10/16/85  11/26/85  112.2 

3  11/08/  85  12/18/85  111.0 

4  12/09/85  01/16/  86  85.3 
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Table  2.      Error   term  calculation. 


Source 
All  plaats 
Rua 
Cross 

Each  plant 


Sum    square  residual 

SSRes(ALL) 

SSResCRUN) 

SSRes(CROSS) 

SSRes(PLANT) 


crtn-m 
r(  c  tn-m) 
c( rtn-m) 
cr(  tn-m) 


^   c   =  No.    of   genotypes,    including   parents;    r    =   No.  of 

runs;    t   =  No.    of   time   intervals   te  s  t  ed/ p  1  ant ;    n  =  No.  of 

test  tubes/time  interval;  and  m  =  No.  of  parameters  in 
the   sigmoidal  model. 


SS(CROSSXRnN)    =   SS  ERROR 

=   SSResC  CROSS)  +  SSRe8(RDN)-SSRe8(ALL)-SSRes(PLANT) 
•^•^•ERROR'      f-SSResCCROSS)  +  ^- ^  •  S  S  Re  s  (  RU  N  ) 

-d.  f.  sSRe8(ALL)    "   ^- ^- S  SRes(PLANT) 
MS   ERROR  =  MS(CROSS  X   RUN)    =    SS    ERROR  /    d.  f gj^j^QR 
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Table  3.  The 

ki 1 1 ing 

time  (T 

ime50)   and  its 

coefficient 

of 

the   standard   deviation  (K) 

for  the 

parents 

and  their 

Fl 

Liy  w  ^  xuo* 

P2 

P3 

1 

76. sy 

62.  3* 

2 

100.8 
6  4.  2 

84.  5 
49.  5 

^3 

111.3 
5  5.4 

114.1 
93.9 

108.3 
49.  5 

113.1 
57.4 

128.8 
54.  7 

81.3 
48.  4 

97.7 
48.  4 

^5 

101.1 
54.9 

140.9 
61.1 

124.0 
54.7 

107.1 
52.7 

123.3 
5  5.6 

^6 

105.6 
50.6 

110.8 
57.5 

117.0 
6  5.5 

106.6 
54.6 

103.6 
58.4 

96.0 
48.3 

^   ?!=   PI   271998,   P2=   PI   324607,  P3=   Oregon  1604, 

P4=  ND  364,   P5=   Wyoming   166,    and  Valley, 
y   Time50,    mean  value   for  4  runs;   LSDq  q5  =  2  8.6  5. 
*    K,   mean  value   for   4  runs. 


Table  4.    Genetic  variance  components   for  killing   time.  < 

Genetic  component  ^  Estimate  and  s  tandard  errory 

D  7  4.  9  +  12  5.5 

Hi  46  8.  7+318.  5* 

H2  435.2+284.  5* 

F  61.9+306.5 

226.0+1  91.5 
B  227.6+47.4** 
2  D  =  additive  effects   of  genes; 
H2^=  dominance   effects   of  genes; 

H2  =  <iominance    indicated    by    asymmetry    of    positive  and 
negative  effects   of  genes; 

F  =  covariance  of  dominance  and  additive  effects; 
2 

h   =square  of  the  dominance  effects  over  all   loci  in 
heterozygous   phase   in  all   crosses;  and 

E   =   environmental  error. 

y*.    **  significant  at   20%,    1%   level,  respectively. 


Table   3.     Hayman's  heritab 

Inhe r i t anc e  parameters^ 
(Hi/D)l/2 

H2/4H1 

Narrow  sense  heritability 
Broad  sense  heritability 

^  (Hj/D)^/2  =  average   degree   of  dominance; 
^2^^^1  ~  average  frequency  of  positive  vs  negative 

alleles ; 

Kq/Kr  =   ((4DHi)l/2+p)/((4DH^)l/2_F)    ^he    ratio  of 

dominant   to  recessive  alleles; 
average  number  of  genes  showing  dominanance; 
Narrow  sense  heritability  =  ( 1 / 4)D/ [  ( 1  / 4)  ( D -F +  Hj ) +  E  ]  ; 
Broad   sense   h  e  r  i  t  a  b  i  1  i  t  y  =(  1  /  4)  ( D-F +  H,  )  /  [  ( 1  /  4)  ( D-F +H,  ) +E  ] . 
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ility  parameters. 

Est  imate 
2.  50 
0.23 
1.40 
0.  52 
0.054 
0.  3  46 


-5- 
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Table  6.  ANOVA  table  of  mean  killing  times  for  the  half 
diallel  cross. 


Source 

Cross 
Run 
Error 
Total 


d.  f. 

20 
3 
60 
83 


Mean  square 
** 


953.76 
4848.77 
410.29 


** 


significant   at   1%  level. 


Figure  1.  Graph  of  variance  vs.  covariance 
6  bean  cultivars  for  heat  tolerance. 
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APPENDIX  A  i 


I 


•  « 

•  Pir  THE  Sltf^CIDAL  fOUtL  TG  ALL  PLANT,  TO  rtbN.  « 
«  rC  CRCSS,   ANO   TO   EACH  PLANT  « 

•  « 

ff 

DATA  all: 

INPUT  PLANT  l-i  TIME  5-7  Al  10-13  I  A2 
A3  20-^3  I  A*  i5-Za  I  A5  30-33  1  Ao  J5- 
RUN  '*0  PI  <i5  P2  50; 
IF  Pl=3  CR  P2=3  THEN  DELETE; 
IF  Pl=l  ANC  f2=i  CR  Pl=l  ANC  P2=l  THEN 
IF  Pl=l  ANO  F2=2  OH  Pl=^  AND  P2=l  fhtN 
IF   Pl=l   ANC    Fi=3   OR    Pl=3   ANO   P2=i  THEN 


15-Id  1 
■38  L 


IF   Pi=l   ANO  P2=4  GR   Pl=4   ANO   P2=l  THEN 

IF    Pl=l   ANO  P2  =  5  QR   Pl  =  !3   ANC   P2=:l  THEN 

IF   Pl=l   ANO  P2=6 
IF  Pl=2 
IF  PI=2 


CR   Pl=6   ANG   P2=l  THEN 


AND   P2=3  CR   Pl=3   AND  P2=2 
ANC    P2-4   QR   Pl  =  i,   AND   P2  =  2 
IF  Pl=2  ANC   P2-5  GR   Pl=5   ANC  P2=2 
IF  Pl=2  ANO   P2=6   GR   Pl-o  ANG   P2=2  THEN 

THEN 
THEN 


THEN 
THEN 
THEN 


IF  Pl=3  ANC  P2-3  GR   PI=3  ANC  P2=3 

IF  Pl=2  ANC  P2=2  QR  Pl=2  ANO  P2=2 

IF  Pl-3   ANC  P2  =  <i   GR   Pl^^r  ANC  P2=3 

IF  PI=3  ANG  P2=5  CR   Pl=5  ANG  P2-J 


IHEN 
THEN 


IF   Pl=3  ANC   P2=6  GR   Pl=6  ANC   P2=3  THEN 


IF  Pl  =  <, 
IF  Pl  =  <, 

IF   Pl=4   ANO  P2=6 


ANC   P2  =  'i   QR    Pl-4    ANG   P2  =  '» 
ANG   P2  =  5   GR   Pl  =  5   ANO  P2='f 
QR   Pl=6  ANG  P2=4 


THEN 
IHEN 
IhtN 

ANO  P2=5  THEN 
ANO  P2=5  IHEN 


CRCSS 

CRCoi 
CRCiS 
CRCSi 
CRCSS 
CRCSS 
CRC5S- 
CKCiS- 
CHCii- 
CRCi;j- 
CRC3S- 
LRCSj^ 
CRCSi^ 
CRCSi^ 

CRCSi- 
CRCiS- 
CRCS5- 
CRCSS^ 
CRC55= 


IF  Pl=5  ANO  P2=5  CR  P1=S 
IF   PL-5  ANC   P2=o  OR  Pi=6 

IF    P1^6   ANC    P2  =  6   GR    P  1  =  6   AND   P2=t>  IHEN 

REP=i;  R=l-(AiyA<i) ;  qgtpui; 

rtEP=2;   R=1-(A2/A5);  CGTPLI; 
REP=3;   R=l-(A3/Abj;  CGTPLT; 
ORGP  Al-Ao; 
CAR03 : 

DATA  cne;  iEi  all: 
PRQC  scrt;  by  plant; 

DATA   ThO;   SET  CNE; 

IF   TIME  >  C   THEN  DELETE; 

PRQC   MEANS   NCPRIM;    BY      PLANT:    WAR  H; 
OUTPUT  OUT=N£h  MEAN=RC; 

DATA  three;  set  nek; 
PRQC  scrt;  by  plant; 

DATA   CCNO;   HER6E  CNE   THREE:    dY  PLANTJ 
IF   TIME  =   0   THEN  DELETE: 
REAO^  l-(R/RC  )  ; 

PRQC  nlin; 

FARMS    B=.Ol    IC    .3   BY    .1   U=50   TG    180   BY  30; 
L  =  EXP{-e«C  riM£-(j  )  )  ; 
MCDEL   «EAC=l/( 1>L) : 
0£R.U  =  -L«B/(  I-fL  J»«2; 
OER.e=<TIME-L)<L/{ 1*L)««2: 
OUTPUT   0UT=MCCALL  R=R£ACALL  PARMS=BALL 


=*1X1« 
-•1X2' 
=•1X3' 
=• IX^' 
=•1X5* 
=•1X6' 
:'2X3' 
=  '2X'«' 
:'2X5' 
='2Xo' 
:'3X3' 
•2X2' 
•iX't' 
•3X3' 
•3Xo" 
"tX^t* 
•tX5' 
•tXb  • 
•5X5' 
•5X6* 
•6X«j' 


GALL  ESS=SSR£SALL 


HKCC   NLIN    ;    EY  Ci*CSS; 

PAHhS   a=.01    IC   .3  UY    .1  U=50   TO   180   dV  30; 

L-Exp(-a«(  iifE-u)i : 

MODEL  READ=l/(  1  +  L) i 
UEK.a=(IIM£-L»«L/(  1+L)0»2; 

CUrPUT  OUT=N£UC  P=PR£AOC  PARMS=aC  UC  tSS=SSR£SC; 
PROC  MEANS  N  NCPRINT:   BY  CROSSiVAR  READ; 

ourpui  ouT=N£bc  n=m:; 

DATA  NEhCJSET   NEhC;   BY  CKCSSi    IF  FIRST. CRCSSJ 
DATA  CKGSSi;   MERGE  NEhC   NEWO;   BY  CKOSS: 
DFC  =  NC-2  : 

DATA  CRaSS2;SEI  CRCSSi; 
PROC   SCRT;bY   CESCENOINto  UC: 
PROC  MEANS  StjM;VAR   SSRESC  OFC; 
OUTPUT   CUT=N£WC   SUM=SSSRtSC  SOFC; 
PHOC   SORT  0AIA=CCNO;   BY  RUN; 

PROC  ^LIN:aY  run; 

P4RMS   a=.01    IC   .3  BY   .1   0  =  50   TO   IBO  aY   30 ; 
L  =  EXP(-B«I IIME-U)  I  ; 
MCUEL  RtAO=l/(l-»Lj; 
D£R.c=(riM£-l,)«L/(  1*LI««2; 
D£R.U=-L»fl/  t  1-»L)««2; 

OUTPUT  CUT=NEWfi   P=PR£AOR   PARMS=   BR  UR  ESS=SSRtiR 
PROC   MEANS  N   NCPRINT;SY   «UN;vaR  RtAO; 
CUTPUI   OUT=N£hS  N=NR; 

UAIA   NthR;SET   N£ki«;BY  RUN;    IF  FIRST. RUN; 

DATA  runi;m£bg£  NtwR  nehs;by  RUN; 

DFR=NR-2  ; 

DATA   RUN2;S£T  RUNi; 

PROC  scri;by  cescenoim,  uh; 

HRCC   MEANS   SUM;VA«   SSKESH  OFH; 

PRCC  print; 

PROC  NLiN;tY  plant; 

PRAMS   H=.Q1    TC   .1   BY    .02   U=   30   TO   150  BY  ^0; 
L^EXP(-d«(  IIME-U); 
MODEL   REAO=l/t UL) ; 
0ER.8=(TIMe-U)«L/l  l+L)»<»2: 
D£R.U=-L«e/{ UL )««2; 

CUTPUI   QUT  =  N£ViA  P  =  PR£AO  PARMS=8P  UP   tSS  =  S5kESP; 

PROC  print; 
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«««««««««««««««««««4««404  WEIGHTtO   ClALLEL  AhALY^iXS  e«<'««»«««««0««w«««<><>«««>««««« 


♦  PAReNT   I  =  PI  27199(J 

« 


2  =  PI   32^160  7  >  .  .0 

3  =  GRfcGGN  ibCr  ^  • 

•  4   =  ND  364  ;•  '  .                                                '        '  «c 

•  5   -  HVCMING    166  :  ^-i  ^ 

•  6  -   VALLuY  « 
«»««04«0«4444444404004444 444, 4„^,^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 

9 

DATA  test; 

INPLiT   RUN   5   TIf£50   10-13   1   STOER   15-20  <,   fSRtS  Z5-iH   5  PI  35  40; 
SIcU=SIOeR/StRI(MSR£Sj;  .      .  ,  V  ' 

SIEU=SI£U«STELi:  -     •  »    ^  ,1 

SIGCP=.06o075: 

cards; 

PROC  scrt;  ey  pi  pz;  ' 

PROC  MEANS  mean;    dY  PI   P2;VAR   TIME50  STEu; 
OUTPUT   0UT=0AIA   MtAN=TIM£50  Ki 
DATA   CATa:    Sfcl   lata;  K=SgRT(K); 

PROC   MEANS  MEAN  CArA=r£ST  NOPKINT;aY   Pi;    VAR   TIMfciU;  '  ^ 

OUTPUT   QUT=NUMe  M£AN=TIME50; 

PROC   MEANS  N  CATA=NUMB  NOPRINT;   VAR  TIMESC; 

OUTPUT   OUT=NUMBER  N=NP; 

DATA  number;   SET  NUMBER;  K££P  NP; 

DATA  CG:  INPUT  L  i  i«;cARi:s; 

CCRR         STE_CORR         PR_CCR=0  ,  '   ; ' 

DATA  a;  input  A  s  n;  cards;  •    -  * 

VARIANCE  CCVAR  h-»V   W-V   Y   YR  STNCRD_Y  PARAEOLA  RANK  CHfl 

DATA  e;  INPUT  e  s  a:  cards;  ~  ; 

0  F  hi  H2   SHCR  £ 

DATA  E;  INPUT  £  S  iii  CARDS;  - 
0_OF_DCM  POSiNEG  OCM:R£C  SENtS  - 
DATA  G;    input   6    i   i«;  CARDS; 

VCLQ   VILI   i,CLCI   VOLI   MLI_MLO  MLI_MLC2   0  F  HI  H2  HSwK  ERRCR 

DATA  x;  INPUT  X  %  ..;  cards;  <: 

INTERCPT  SLCFEtej  B_£RROR 

DATA   go;    INPUT   RUN    S   a«;CARDS;  "        .  ,) 

l«-V_hCM        PR_h-V_HCM        b-C        PR_a  =  0        B=l      PR   B=I  •-' 
DATA  ERR0R;S£T   TEST;  . 

IF  _N_  =1;  k 
KEEP  SIGCP; 

PRCC  matrix; 

FETCH   P  DATA=:NUMBER;«««*CATA   set   klTH  NUMBER  CF  PARENTS  GNLY: 

IVrru   c   rVA"^ctlt'     ^^^""^^TA  SET   i.ITH  PI,   PI,  ObScRVATION  AND  -tlGHIIN*.  VALUtSi 

FETCH   E   DAIA=EfiRCR;   ««0«CATA   SET   KITh  fcRRCR   TERM  GNLY;  ««t.uci, 

U=J(P,P,0):   ♦••♦OBSERVATION  VALUES; 

K=J(P,P,0)  ;«»««K  VALUES  FOR  hcIGHTINS; 

UR=J(P,I,0); 

AR=J(P,1,0); 

SR=J(P,1,0); 

VR=J(P,1,0) ;«««*VARIANCE    IN   AN   ARRAY;  '  ■* 

hR=J(P,l,0);««*«CQVARlANC£    IN  AN  ARRAY; 
•BUILD  OlALLEL  DATA  ARRAY; 
U=NRGta(C); 

•DETERMINE  NLMEER  OF  MEANS   IN  OlALLEL; 


33 

DO  1=1  IQ  xk: 

P2=UCIf2);  o««op2,   PARENT  <2; 

U«PltP2l=ClIt3J ;  ««««Ut    OilStRVARTIQN   VALUfcS  MA7HIXJ 

uiP2.Pi)=G<i,3);        ««««u,  oasfcRVARriQN  vALuei  matrix; 
KiPi,P2)=c(it<i):       04<>«K.  kiEiuHT  matrix; 
KipatPii=e(it^j;       ««««Kt  weight  matrix; 
end: 

print  g  u  k; 

«CIALLfcL   ARRAY   IS  MATRIX  U: 

UH  =  VeCOIAG(U)  :  parent  OaSERVATIOh  VALLtS; 

KQ0=SSG(KU/(PSP);  AVERAGE   MALUc  Uf  WEIGHT; 

E=E»KOa;  ERROR  CaRRfcCTEU  FUR  MtAN  VALUE:*; 

MI  =  VEC0IA6(K)  ;  ««««  PARENT   xEIGHT  VALUtS; 

UP=SUM(UIIS/ (KIXsKII ) iU/SSQi  iU/Kll ) : 
PRINT   UII   KQG   t   Kll  UP; 
OC  1=1  TO  p; 

UK( I fl  )  =  SUM(L( .1 )S/(K( «I )aKt ,1 ) ) )tf/SUM( (  1 4/ t K ( , U ( . I ) )  )  1  ; 
ARI  I,I)=SUM(li(  1 1  )S/  (K(  ,1  )SK  (  ,1  )  )  1  ; 
SRI  I .1  l  =  SUM( Is/ (K( «I ) tK( T I)  ) ) : 

VR(  I  ,n  =  (  IS/ (p-n  isissui  luc ,  iis/ki «i i  > i-i lri  i,i)»uki  it 1 1  usumi  (  is/(k( ,i  ).»k,i .ii  i 

)  DSKOC; 

l«RU  f  I  )  =  l  la/  (P-U  )»(SUM<  I  IUII-UPI*/KII)S(L(  f  I)-URi  1,1  J  )J«/K(,i)  )  isKOC; 
ENO; 

PRINT   UR  AR   5R   VR  WR; 

VOl.C=(  IS/ t  P-  U  mSUMl  HjllB/KIIJSIUIIa/KII  ))-(UP#UPsSU«t  la/  (KI  IdKIin  I  )-K00; 

viLl=SLiMC  VR(  ,1 )  »a/p; 

WCLQI=SUM( HR ( , n ]S/P: 

VOL  1=1  la/IP-n  1S13SGIAR)-<SUM(  AKi<SR)««»2»/SSulSRJ  )  J  BKCO-KOC«/ I  PaP  »  ; 

MLI_MLC=t  <SUM<Ua/(KSKnB/SUMI  IS/I  KaK)  )  l-(  SUM  I  U 1 1  J/ ( K  1 1 JK  U  )  la/SUM  I  la/ ( K  I  UK  1  i  )  )  1 

) ; 

Ml.I_ML02=MLI_MLC««2; 

n=p: 

print  vqlo  vili  hclqi  voli  ml1_mlq2  mli_mlc  n  pi 

oc«=vcLO-6;  ~  ^:  ' 

F=2*VGL0-4«taCLCI-l (2«(N-2 J l*E)a/N; 

Hl  =  VCLC-<i<>hOLCI  +  4«VILI-(  <  1 3«N )  - 2  J « 1 1  J/N  ; 

H2  =  't«VILI-'»«VCLI-2«t; 

MSiJR  =  4«MLI_MLC2-('»«iN-l)«£)  s/N««2; 

DEG_DCM=St;RT  Ihla/bUMl  ; 

P0S_NEG=H2S/(<t«Hll  ; 

DGM_R£C=CS4,RIl4»CQMOHU+f  la/JSuRTC^UCMOMIl-Fl; 
G£NES=HS£Ra/h2  ; 
CDaM=(N»«5+N««4lB/lN««5l  ; 

Cf  =  (  ('.«(N«*5  )  )  +  <20«  (N««^) )-(  16«(N««3)  l-f  i  1 6« (  N*<»2 )  »  I  a/<N««5)  : 
CHl=(N««3'f  ('>l«(N««<r)  l-(12«(N««3)  )«<4«IN««2n  )S/IN««5i; 
Ch2=  «  Sb^lN***!  J  IS/(NO«5l  ; 

CHSaR=(  {  l6«(N*«'il  l  +  (  160(N««21  )-(32«N)->16la/(N««5i  : 
C£=IN««4)»/(h««5); 

n_v=hR-vfi; 

»<_PLUS_V-=hR+VR  ; 

«  STANCARO   EBBCRS  i 

S2=(  iSSCiIta_V)-((SUM(K_V))»»2)»/N)a/HN-l)82)  i: 

SCOM=SGRT(  CCCCJ">S2  J  I  ; 

SF=SCRI(  ICF«S2H  ; 

Shl=SCHT ( ICH1«S2 J )  ; 

Sh2  =  SQRTI (CH2«S2))  ; 

ShSQR  =  SQRTUChSGR«S2J  i: 

SE=SGRT ( (C£«S2 ) ) ; 

EBROR_CF=   SCCMI  I SF I  I Sh U 1 SH2 I ISHSGkl ISE; 


o  T   TEST   SIAST   FOR  0   F  HI  H2  HSu  ERRUR;  ^ 
DCNTT=CQ«»/SCCH;  ^ 
IF  OCMTT  >  a    THEN  OT=-l«tCMTT; 
IF   DCMIT  <   0    IhEN  0T  =  C0f1Ti; 

ftt=f«/sf; 

IF  FTT  >   C  ThEM  FT=-l«FTT; 
IF  FTT   <  0   THEN  FT=FTT; 
HiTT=Hl»/SHl ; 

IF  HITT  >  0   THEN  H1T=-1*HITT; 
IF  HUT  <  C   THEN  H1T=HITT; 
HiTT=H2a/SH2 ; 

IF  H2IT  >  0   THEN  H2T=- l'>H2 TT  J 
IF  H2TT   <  0   THEN  H2T=H2TT; 
HSaRTT=HSQR«/ShSGR: 

IF  HSQRTT  >   C   THEN  HSCR T=- l«HSURTT ; 
IF  HSQHTT  <   C   IhEN  HSGR I=hSORTI ; 
£IT=£3/S£; 

IF   ETT  >  0  THEN  £T=-i«ETT; 
IF   ETT  <   0   THEN  £T=£TT; 
PCT=P«CBTICT,N-l); 

po=2«pct; 

PFT=PR08TIFT,N-l); 

pf=2»pft; 

PhlT  =  PRQeTIHlT»N-l  J : 

Phl=2«PHlT;  \ 

PH2T=PR0fiT(H2T,N-I j; 

PH2=2«PH2l; 

PHS(3RT  =  PR08T  IHSCRT»N-I)  ; 
PhSQR=2«PHSCKT ; 

PET=PRCaT(£T,N-lj;  ' 
PE^2«P£T; 

PRQa_T=P01 IPFl IPHil 1PH2I 1 FHSGRI IPE; 

ncte  oiallel  cross  data; 

PARfcM  =  UII'; 

SIOCCV=«SSa(hR)-<  (SU«lkiRJ«*2)«/NI  )«/lN-lJ  ;     K  ■ 
STQPAR=SCRT(VCLOaKOO) ;  f    .  ^' 

ST0VAK=<SSC1  VR)-I  (SUMIVR  J«'»2)«/N)  )B/(N-IJ  ;  '  * 
PM  =  SlJf(<PAR£NI  )«/N; 

Mn=sun(MR)s/N; 

v«=su«(VR  j«/n;  .,  -  ■ 

YR=uii'«/Kii';  '  f   :    '  ^ 

SI0_Y=(PAR£NT-CP)»/1SIDPAR»KXI');  -       .     :  5 

KO=RANK(h_PLLS_Vi; 

PARAbOLA  =  Si;RT(  VRaVQLQ): 

VR1=WR'  ;kRl=t,R'  :i,_PLUSV  =  h_PLOS_V  ; h_V !=>._  V  ; PAKAocL  1  =  PAR AfaCLA '  ;   RCl  =  RO'  : 
STATS  =  VRl//i,Rl//W_PLUSV//W_Vl//PAREM//YR//STD_Y//PA«AbCLl//KGi; 
FETCH  C  OArA=A  TYPE^CHAR; 

NOTE  ARRAY  STATISTICS;   PRINT  STATS  RCliNA«E=C: 
PLOT=STATS'; 

OUTPUT  PLOT   CLI  =  PUT1    ( RENAME=  ( C01.1  =  V AR I ANC£  C0L2=C0VAR  CCL3=I.  PLUS  W 
COL^=H_V  CCL5=Y  C0L6=YR  ~  ~ 

CCL7=STNCR0_Y  CCL8=PARAeCLA  CaL9=RANK  ORO  ROH=PLr H . 
CUTPUT   STATS   CliT  =  PLT; 

MEANS=CEG_OCM I PCS_N£G N DCM_R£C I  I  GENES; 
NCTE   HAYMANS   ANALYSIS  QUANTIESS 

0UANTITY=VOLCI  I  V  IL II  1  WQLC 1 1  I  VOL 1 1  I « L I_MLQ J 1 MLI.MLU^  I  lOOHl  J F  1 IHI I  |h2  I IHSCRI l£; 
CUTPUT  QUANTITY  CUT=VCLC   ( R£NAM£= ( CCL 1= VO LO  CCL2=VILI   CUL3=I.ClOI  CaL4"vCLI 
CCL5=«LI_HLC  CCL6=NLI_MLC2  CCL7=D  CQLa  =  F   CaL9-Hl  CCLI0=H2   CCLll  =  hSi.R  Ct)Ll2=ERRUR 
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FETCH  Y  OATA=G  rYPe=CHAft; 


PRINT   tUANTIIY  CCLNA«6=Y; 
FfcTCH  X7   CATA=e  IYPe=ChAR; 

NCre  STANOARC   errors;   print   tRRUR_Qf  CULNAMt=X7; 
NCIE   T  PRa£=C;pRINT   P«UB_T  CCLNAMe=X7; 
FETCH   I  CAIA=E  IYP£=CHAK; 

NOTE  MEAN  EFFECTS  OVER  ALL  PARENTS;  PRINT  MEANS  COLNA»(E= 
MhV  =  SUM<  Ii_PLLSV)«/N; 

Cl=SUM I ( H_PLLSV-M«W la (PARENT-UP JJ/KI I • ) ; 
C2  =  SSQ(W_PLUSV  )-( IS0MIW_PLUSV)*»2)«/N) ; 

C3=vOLOsK00a (P-1 ) ; 

CCRR=C1»/SGRTIC2«C3J ; 

SIEC  =  SC,RT(  <  l-CCRR««2}*/«h-2J  J  ; 

CCRT=CCRRS/ST£C: 

IF  CCRT  >   0   THEN  CT=-I«CCRT; 

IF  CCRT  <  0   THEN  CT=CCRr; 

PRG8CA=PRC£I (CT ,N-2) ; 

p«0BC=2«p«cacA : 

CCRR1  =  CCRR1 ISTECI  IPRCfiC; 
FETCH  CO  OAIA=CC  TYPE=CHAR; 

NOTE  CORRELATION  OF   Y  ANO   H+V;   PRINT  CCRKi  COLNA«t=CC; 
VV_Vki=(STa\(AB-STCCCV)««2; 

vvxvh={siovar«stccqv)  ; 
pcov^vrshR; 

CVRKR2  =  I  «SU«(PCCV)-{  (SUMIVR)«SUm  hR)  )S/M  la/ IN- 1)  ; 
T2=l  (N-2)»VV_\(li)«/(  t  VVXVh-CVRKK2)«4)  ; 
Bl=SljmVH«hR)-«  (SUM<  VR)«SljMChR)  )a/N)  ; 

B2=ssaivR-y(M) ; 
a=bis/B2; 

Se  =  S(.RT{  ( (STCCCV«(N-U  )-(a«BlJ  la/CIN-2)»B2J  »  i 

oiF_o=e»/S6; 

0IF_1=( L-8»a/S£; 

IF   12  >  0  THEN  I2A=-1»T2; 

IF   T2  <   0   THEN  T2A=T2; 

IF  OIF_0  >  0    THEN  0 IFC=- 1«0 I F_0  ; 

IF  DIF_1   >   C   THEN  0 1 F 1=- 1«Q I F^I  ; 

IF  OIF_Q   <  0   THEN  OIFC=OiF  0;~ 

IF  0IF_1   <  0    TFEN  OIFl=OIF_i: 

PR0aT2A=PRCBT (  I2A»N-l); 

PHQaT2=2«PR0£T2A; 

PR06AC  =  PR0£T (CIFCtN-2  1 ; 

PROaC  =  2»PROBA0  ; 

PROBAl  =  PROeT (CIFI,N-2)  ; 

PRaBl=2«PRCBAi; 

INT£RC£P=hf-£«VH; 

NOTE  GRAPH  STATISTICS; 

GRAPH=INItRCEP  1  I  El  I  SB ; 

FETCH  XI   OATA=X  TYP£=CHAk; 

OUTPUT  GRAPH  CljT  =  LINEl  «  R ENA«e=  « COH=  INTE BCP T  C0L2^SL0PE 
PRINT   GRAPH  CCLNAHe=Xi: 

TTEST  =  T2I  IPRCET2J  IDIF_0|  I  PROEO I  I  0 1 F_ 1 1  I PRC£ I ; 
FETCH  X2   OATA=GC  TYPE=CHAR; 

NOTE   T   TEST   STATISTICS  FOR;    PRINT   TTEST  CCLNAMc=X2: 
NOTE  NARRCh  SENSE  HERITABILITY  ESTIMATES; 
HRITELTY=JQ.25«C0«14/»  (0.25*«0QM-»Hl-F»  l+c  Ji 

PRINT  hritelty; 

PHOC   CORR  DATA=PLTi; 
VAR   YR  h_PLUS_v; 
PROC  REG  DAlA=PLTi; 
MCOEL  CQVAS=VARIANC£; 
DATA  PLTi;S£T  FLTi; 


IF  PLI  =  'RQI«1'  IhEN  OO;  C  l^CQ  W  AR  ;P  LAN  f  =  •  P  I  27l99d  'itUbi 

If  PLT  =  'RCh2'  THfcN  OG;  C2  =CC  V  AR  ;  P  LAN  r=  •  P I   J2'i607  •;£MJi 

If  PLT  =  «RCW3'  IHEN   DO:   C3=CG V AH  ;P LAN  I- • OREbON  1604';tNUi 

IF  PLT='RCl.^'  IHEN  DO;  C4=C0 V AR ;P LANT=  •  NO  3o^  'JEMiJ 

IF  PLT  =  'R0ta5*  IhEN  00;  Cb=LO  VAR  ;PLAM=  •  k  Y  CM  I NG  ibb';£NU; 

If  PLI  =  'RQ1.6'  THEN  DO;  C6=C C V AR ;P L ANT  =  •  VA LLE Y  •;ENUi 

DATA  parab;  set  volo: 

Xl=0-5+(0.5«SGRT  (l-I  i  ^«  » kiOLD  1- V I L 1 )  l/VQLQ  J)  I  ; 
X2=0.5-<0.5«SCRTU-ll4«thOLOI-VILm/VCI.O  ll»  J 
V01=VCL0«XL««2 ; 

wci=vcLO»xi; 

VRl=V0L0«X2««2;  I 
WR1=V0L0«X2;  ' 
INrER=llQLQI-VIH; 
kiPV  =  hOLOI*VILl; 

PRQC  print; 

DATA   PARA8i;S£T  PARAB; 

x=voi;yi=woi;var=o;ccv=inier;ootpui ; 
x=vRi  ;Yi=xfii  ;vAR=wRi  ;cDv=*iRi;aurpui  i 
x=vili:yi=iiClci;var=.;ccv=.;  quipui;  - 

KEEP  X  Yl   VAR  COV; 

DATA  RUNi;  SET   VCLO:  ^  .. 

00  VARIANC£=0   TC  500   BY  5; 
PAR82=-i«(SCRT< VARIANCE«VCLO» ) ; 

ooipijT  ;£No; 

KEEP   PARB2  VARIANCE; 

PRQC  SCRT;   fiY   CESCENOING  VARIANCES 

DATA  L2;S£I   LINEi;  ,  I 

00   VARIANC£=Q.CO   TO  500   6Y   5:  " 

L2=INTERCPT   ♦   « SLOP£«VARIANCE  )  i 

OUTPUT  :£N0 ; 

DATA  R0N;S£T  VCLC; 

00  VAfilANC£=O.CO   TO   500  flY  5i  ]^ 
PARa2=SQRTIVARIANCE*VCLQ); 

OUTPUT  :£No;  -V    ;  ' 

KEEP  PARe2  variance; 
PROC  sort;  by  variance; 

DATA  RUN2;S£T   RUNI  RUN; 
DATA   PLCT;S£T   PLTl  RUN2: 

LAb£L  ci=ccvariance;  '"  • 

PROC  REG  CUTES7=£ST  0ATA=PLT1; 
MODEL  YR-W_PLUS_V; 

PROC   MEANS  MEAN   OA rA=PLT I ; V AR  YR; 
OUTPUT  QUT=M£ANY  MEAN=YR; 

DATA  EST;  SET  EST ;b=w_plus_v ;  i  •     -      .  •■'  ■ 

t£st=i:k££p  TEST  B;  jr  \  ■     »'  i 

DATA  PARUSET   PARAB:  <  ?  1       i  J 

KEEP   UPV ;  .  .  ,  / 

DATA   LINE;MERG£   EST  MEANY  PARI;  *"      '  '      ^  ' 

PROC  sort:  by  test: 

DATA   t<R;S£T   PLTi:TEST  =  i:  - 

yri=yr: 

keep  k_plus_v  test  yri: 

DATA   MRi:   SET  FARA8; 

w_PLUS_v=taoi*vci:rEST=i:cuTPuT;  -  » 

w_plus_v=i.ri-»vri:test=i;cuiput:  .  ' 

KEEP  test  ii_plus  v:  • 

DATA   kR2;SEl   kR  hRi: 

PROC  SORT:   BY  TEST; 

DATA  GC:M£R6E   IiR2  LINE   :aY  TEST; 
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6SI  =  yR+«d«(ki_FLUS_\(-wPV)  1  i 

PRGC  print; 

PRGC   print;   VAR   YRl    EST   YR   fc   W_PLlJS   V  WPV; 
OAIA   GPLCT:    set   PLCT   PARAfci  L^; 

KctP   VARIANCE   PARa2  A   Yl   VAR  COV  CI   C2   C3  Ct  Ca   Co   PLANT  LZ; 
PRGC  tiPLCT: 

TITLEl  VR/hR  CIALLEL   6  PARENTS; 
SYM8CL1   C=BLACK   V=l  F=SI«PLcX; 
SYMaCL2   C=BLACK   V=2  P=SIHPLex; 

SYMBCLJ  C  =  bLACK   V=3   F  =  SIMPLEX;  ,.  ^  ■»  t  ,  ,^  . 

SYMBOL*  C  =  eLACK   V=4   F  =  Sl>'PLtx;  ■'  .    .  [  r  ^ 

SYMBOLS   C=eLACK   V=5    F=SIMPLtX;  ;  . 

SYMB0L6  C=eLACK   V=6  P=SIMFL£X; 

SYMBOLB  C  =  aLACK   V  =  M  H=2   F  =  SPeCIAL  L=L   I=JCIN;  '  ■* 

SYMBCL9   C=aLACK   L=l  I^JOIN; 
SYMBCLIO  C=aLACK  L=I  I=SPLIN£; 
SYM8CLU   C  =  ELACK  L  =  3    I  =  JCIN; 

PLOT   CI»VARIANC£=l  PARB^OVAR IANC£= 10  C6»V ARI ANC£=o 

C20VARIANC£=2   C3»V  AR I  ANC£  =  3   COVARIANCE^*  l.5«VAKIANCc  =  3   L2«VAK1  ANCc  =  9/0 


«««««««« i>««««««>)«o«««««««o«««i:i««««««<««««««9«o««<i«i>o« 
e  ««««««««««««««««««  ^       ^  4^^^^^^  ^j,^^^*^,^^^^^^^^^^^^^^ 

«  « 

♦  DATA  FOR  OIALLEL  ANALYSIS  • 

oo« «««««««« o«*c«««««o«»«««*««««««i^««««^«*««^«^^e^^^^ 4 
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Abstract 

To  estimate  the  genetic  component  of  cellular 
membrane  thermostability  in  Phaseo lus  vulgaris .  parental 
and  Fj^  plants  from  a  6-parent  half  diallel  were  tested 
by  electric  conductivity  and  the  results  were  analyzed 
by  Hayman's  method.  Membrane  thermostability  was  found 
to  be  a  quantitative  trait,  with  environmental  and 
dominant  effects  accounting  for  most  of  the  phenotypic 
expression.  Narrow  sense  heritability  was  relatively 
low  (5.4Z) 


